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^ , We explore the thermodynamics of the gluon plasma in SU(3) Yang-Mills the- 



> 

o 

(N 



ory emerging from the non-trivial spatial dynamics of valence gluons. The lattice 
data suggest that these gluons interact with each other linearly at large spatial sep- 
arations. At high temperatures, valence gluons should reproduce the pressure of the 
non-interacting Stefan-Boltzmann plasma along with the leading perturbative correc- 



jil^' tion. These properties of valence gluons can be modeled in terms of the integral over 

1-^. their trajectories. We calculate such a world-line integral analytically and obtain 

^ \ the pressure and the interaction measure (e — 3p)/T^ of the gluon plasma. Addi- 

> ■ 

tionally, we account for the contributions of stochastic background fields to these 

O " 

CN ■ thermodynamic quantities. The results turn out to be in a good agreement with the 

(N 

QQ . corresponding lattice data. In particular, the lattice-simulated peak of the interaction 

o 

00 ' measure near the deconfinement critical temperature is reproduced. 

o ■ 

[x^ I. INTRODUCTION 

M'. 

How strong are the interactions in the quark-gluon plasma? In connection with RHIC 

n 

heavy-ion experiments, this question has been posed p| to understand the low shear viscosity 
necessary for a hydrodynamic description of the early phase of a heavy-ion collision |[2]. On 
the lattice, strong interactions in the QCD plasma up to temperatures about twice the decon- 
finement critical temperature are known from simulations of the interaction measure ^, ^. 
Analytic models aiming at a description of these lattice data use both perturbative [5| and 
nonperturbative jsl techniques. In the present paper, we revisit this problem by analytically 
calculating the pressure and the interaction measure {e — 3p)/T^ of the gluon plasma. The 
main nonperturbative dynamical input of our calculation is the spatial confinement of va- 



lence gluons and the Polyakov loop, w 
respectively. Following previous works 



lich are known from lattice simulations 



7| and 



^, we use the term "valence gluons" emphasizing 
their confining interactions at large spatial distances. At small distances, valence gluons be- 
come quasi-free. Accordingly, the spatial Wilson loop of a valence gluon obeys the area law 
at large distances and goes to unity at small distances. Additionally, we model small spatial 
Wilson loops in such a way as to reproduce the known leading Oi^g^) perturbative contri- 
bution to the pressure of the gluon plasma [l^. We manage to analytically calculate the 
world-line integral representing the free energy of a valence gluon at any distances. Further- 
more, we account for the interaction of spatial valence gluons with temporal background 
fields through the Polyakov loop and for stochastic background fields, which provide 
spatial confinement of valence gluons. It turns out that our model is able to reproduce the 
lattice data on the pressure and the interaction measure of the gluon plasma rather well. 

fin 

The world-line integral representing the one-loop free energy of a valence gluon p, lUj 
can be reduced to the Euler-Heisenberg effective Lagrangian in an auxiliary Abelian field 
of a purely geometric origin. The weight factor associated with this field is constructed 
in such a way as to reproduce the Wilson loop in the original non-Abelian theory. We 
emphasize that we are not considering a YM theory in a constant Abelian-type background 
field jsG]. Another related approach to incorporate nonperturbative contributions to the one- 
loop effective action [l^ uses the two-point correlation function of gluonic field strengths 



provided by the stochastic vacuum model 



The paper is organized as follows. Section II is devoted to the calculation of the world- 
line integral for a valence gluon using a parametrization for the minimal area in terms of the 
contour of the spatial Wilson loop. In Section III, we numerically evaluate the pressure and 
the interaction measure of the gluon plasma. In Section IV, we summarize our approach and 
the main results of the paper. In Appendix A, we evaluate the contributions of stochastic 
background fields to the pressure and the interaction measure. 



II. THE SET-UP FOR THE FREE-ENERGY DENSITY 



Lattice simulations 



in the deconfinement phase (T > Tc) show that averages of 



large spatial Wilson loops still exhibit an area law. In SU(3) YM theory, the critical tem- 



perature is Tc ~ 270 MeV |^ (see e.g. Ref 



17l | for a review). This "magnetic" or "spatial" 



3 



confinement of spatial Wilson-loop averages is the main established nonperturbative phe- 
nomenon in the quark-gluon plasma. It does not contradict the genuine deconfinement of a 
static quark-antiquark pair represented by averages of space-time Wilson loops at T > Tc. 
We follow the strategy of Refs. 0, y, to model the spatial confinement of the valence 
gluons by stochastic background fields B'^. With the ansatz A'^ = B'^ + a'^, the Euclidean 
YM partition function at finite temperature has the following form: 

Z{T) = (^jval{^,t)e^V f^dx, j^Sx{F;M]? )• 

Here /3 = 1/T, V is the three-dimensional volume, -F^^^i^] is the QCD field-strength tensor, 
and the averaging over the background -B^-fields is symbolized by the brackets (...) = {■ ■ 
All the functional integrations imply periodic boundary conditions a^(x, /5) = a^(x, 0), 
-B°(x, /5) = -B°(x, 0). Performing Gaussian integration over the valence gluons, one obtains 
the gluon determinant: 

Z{T) = {{det[-iD;[B]r]y'^-'^''^-'^) , (1) 

with the covariant derivative {D^[B]f^Y = d^f^ + f"''^''B^^f^. Equation ([T]) includes the color 
degrees of freedom of the valence gluons and their physical polarizations, 2(A^^ ~ !)• In the 
one-loop approximation for the a^-field, it can be simplified further: 

Z{T) = (exp {-{N'^ - l)Tr In [-{D^lB])'] }) ^ 

^ exp {-{N'^ - 1) (Tr In [-^^[i?])^])} . (2) 

The loop of fast field fiuctuations receives infinitely many contributions of the slow S^-field. 
In Eq. (l2|), "Tr" includes the trace "tr" over color indices and the functional trace over 
space-time coordinates. 

The free-energy density F(T) is defined by the standard formula 

(3VF{T) = -\nZ{T). (3) 

The interactions of valence gluons with the background fields can be calculated through 
the world-line method. Representing the trace of In [— (D" in terms of the Schwinger 



proper time s as Tr InO = — ^ Tr e one has 19| 
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F(T) = -(iV^ - 1) . 2 V r - / Vz,e-yo^^^l {w[z^]) . 



(4) 



Each world-line path a;^(r) has a center = 7 Jq drx^^r). The integration over the positions 
of the centers yields the factor PV on the left-hand side of Eq. (I3|). The vector-function 
Zfj,{T), T G [0, s], in Eq. (JH) characterizes the shape of the contour of the Wilson loop relative 
to the center c^. The integration goes over world-line paths ^^(r) = x^(r) — obeying the 
periodic boundary conditions 

Zi{s) = Zi{0) + pn, z(s)=z(0). (5) 

Furthermore, the summation goes over the winding number n with a factor of 2 accounting 
for the modes with n < 0. The zero-temperature part of the free-energy density correspond- 
ing to the zeroth winding mode, n = 0, has been subtracted j^, [2^. For simplicity, we use 
the standard approximation which neglects the gluon spin term. Equation ([4]) represents 
the main idea of our approach that the valence gluons, as high-energy fluctuations of the 
gluonic fleld, propagate in a low-energy stochastic background, which enters Eq. ([4]) in the 
form of the Wilson-loop average 



{W[Zf,]) = (tiVexp (^i^ dZf^B^ 



(6) 



We use the notation = B^t^^^ with (tadj)'"^ = —if"'^'^. For brevity, we will call the 
average ([6]) just "the Wilson loop" {s?]. 

We will argue that the Wilson loop factorizes into a purely spatial Wilson loop and an 
averaged ra-th power of a time-like Polyakov loop: 

+00 

{W[z,]) {W[7]) n {L-{T)), (7) 

n=— 00 

with 



and 



= (tr P exp (^^ ^^^^ da,,F,,[B]^ 



{^{T)) = (tiTexp [in dziB^ 



(8) 
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In Eq. ([7]), we have rewritten the spatial Wilson loop through the non-Abelian Stokes' 
theorem. We have suppressed for simplicity the Schwinger strings, which connect the point 
where the field-strength tensor Fj^^ is defined to a given reference point on an arbitrary 
surface S(C) encircled by the contour C. The integration in Eq. ((HI) goes over the purely 
spatial part of that surface characterized by the differential elements dajk- By using the 
cumulant expansion, one can express the Wilson loop in terms of the connected 

averages (cumulants) of chromo-magnetic and chromo-electric fields, = |ejjfcFj\,[_B] and 



El = iF^- [B] . Lattice data [21| suggest that the cumulant-expansion series can be truncated 
at the second term, since the amplitudes of higher terms are significantly smaller. A general 
Wilson-loop calculation leads to three correlation functions 

{E,{x)Ek{x')) , {H,{x)Hk{x')) , and {E,{x)Hk{x')) . (9) 

Following the lattice data [2II, we assume that the correlation function {Ei{x)Ek{x')) van- 
ishes right after the deconfinement phase transition. This fact ensures deconfinement. The 
same lattice data also indicate that the amplitude of the mixed chromo-electric-chromo- 
magnetic correlation function {Ei{x)Hk{x')) is by an order of magnitude smaller than the am- 
plitude of the chromo-magnetic-chromo-magnetic correlation function {Hi{x)Hk{x')) . This 
property leads to factorization, Eq. ([7]). After factorization, the world-line integral over 
Z4{t) becomes a one-dimensional world-line integral for a free particle. Carrying it out, one 
arrives at the following expression (cf. Ref. {^): 

F{T) = -{N',-1).2J2 / -^7=^(^"(T)) /Pze-Ho^^^^iy[z]). (10) 

The dynamics of the strongly interacting gluon plasma is now encoded in the Polyakov 
and spatial Wilson loop expectation values, which make the free energy nontrivial. An 
expectation value of the spatial Wilson loop (W^[z]) depends on the area S = S[z] of the 
minimal surface defined by the spatial contour z(s). At small spatial separations between 
valence gluons, (VF[z]) 1, that corresponds to a non-interacting Stefan-Boltzmann plasma. 
At large separations, a nonperturbative linear interaction between valence gluons emerges 
from the spatial area law, (Vr[z]) exp[— ^(T)^]. By "interaction" we mean the "potential" 
for the gluonic pair evolving in some spatial direction. The linear "potential" between valence 
gluons at large spatial separations can originate from their interaction with chromo-magnetic 
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background fields H^. 

Furthermore, we mimic the leading 0{g'^) perturbative contribution to the pressure [l^ 
by modeling perturbative interactions of valence gluons in terms of the spatial Wilson loop. 
The corresponding short-distance ansatz for the renormalized spatial Wilson loop reads 



e 



C{T)g^{T)TL 



where L = dT\z\ is the length of the contour, and the prefactor C{T)g'^{T)T is chosen in 
such a way that the leading C((?^)-correction to the Stefan-Boltzmann pressure is reproduced 
correctly. We associate the 0{g^)-ierm. in the pressure with perturbative interactions due 
to the massless-gluon exchange at short distances < 0[^^^. These distances are smaller 



than 0{jf) and 



(cf. the recent paper 



perturbative one 



, where the gluon mass terms in the propagator become effective 



22|). First, by comparing our result for the pressure with the known 



10|, we will determine the limiting value C of the function C{T) at T ^ Tc 



This limiting value alone is not sufficient, and in the second step, we will parametrize the 
functional form of C{T). 

Our short-distance ansatz for the spatial Wilson loop can be unified with the spatial-area 
law to a single formula 

{W[z]) = exp \-C{T)g\T)Ty/S - (t{T)s] . (11) 

For the short-distance part, we have estimated the length of the contour, dT\z\, as of 
order Note that, in this way, just one single equation (fTTl) reproduces both the small- 
and the large-distance regimes of the Wilson loop. Indeed, the first term in the exponent, 
which describes perturbative interactions of valence gluons, dominates for sufficiently small 
areas, whereas the second term dominates for sufficiently large areas. 

We implement now an ansatz for the minimal area S, which is crucial for the subsequent 
calculation: 

1 r 

S ^- dT\z X z|. (12) 

2 Jo 

It parametrizes a parasol-shaped surface made of thin segments, as depicted in Fig. [TJ Note 
that, since drz = 0, the point where the segments merge is the origin. Thus, ansatz ffT2!) 
automatically selects from all cone-shaped surfaces having z(r) as a boundary the one of 
the minimal area. We further approximate the integral in Eq. ( fT2ll as follows: 

5 ~ Vf2 (13) 
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Figure 1: A view from above on a typical trajectory in the world-line integral, Eq. (iTOl) . depicted 
by the thick line, and the associated minimal surface parametrized by Eq. (fT2]) . depicted by the 
thin lines. The point where the pieces merge is the origin. 

with 

1 r 

f = - / driz X z). 

The left-hand side of Eq. (fT3l l. given by Eq. (fT2ll . can be larger than its right-hand side. 
This happens if, in the course of its evolution in the spatial directions, the valence gluon 
performs backward and/or non-planar motions. Should this happen, the vector product 
(z X z) changes its direction, and the vector integral on the right-hand side of Eq. (fT3l) 
receives mutually cancelling contributions. The non-backtracking approximation adopted in 
Eq. (fT3ll is widely used in the context of minimal surfaces allowing analytic calculations of 



the associated world-line integrals (see e.g. We further proceed with calculating 

separately the perturbative part and the full free-energy density of the gluon plasma. 



A. The perturbative contribution to the free-energy density 



At T > 3Tc, the Stefan-Boltzmann law along with the leading 0{g^) perturbative correc- 
tion io| . 



45 



15 
IQtt' 



■g\T) + 0{g' 



(14) 
There- 



reproduces more than 95% of the lattice-simulated pressure of the gluon plasma j^. 
fore, at such temperatures, one can retain only the first term in the exponent of Eq. (fTTl ). 
where C{T) should go to a constant C with the increase of T. We determine the value of C 
by calculating the perturbative contribution to the free-energy density with our ansatz for 
a small Wilson loop. 



-C32(T)T(f2)i/4 



(15) 
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and equating the result to the known expression, Eq. ([T4[l . The perturbative contribution 
to the free-energy density has the form 



n=l 



7 _ '-' 

as e 4s 



s Va 



ITS 



Vze-y^dri' {W[z]) 



(16) 



where we have used the leading-order perturbative result for the averaged Polyakov loop 

{L{T))^ = l + 0{g'). (17) 

To calculate the world-line integral in Eq. (fTGl l. we disentangle the fourth root in the exponent 
of Eq. ( fTSl ) by introducing two identical auxiliary integrations as follows: 

{Cg\T)TyP' 



iw\ i\ 1 r r 



—A — 



64A2/i 



We further use the Hubbard-Stratonovich trick 

1 



(47rA)3/2 

where A > 0, and H is an auxiliary constant Abelian magnetic field. It yields 



(18) 



64 



1 



dfifie / (i^ifexp 



16A2/i 



iCg'iT)Ty 



rH^ + im 



7r5/2 {Cg'{T)Ty J, 

This expression, plugged into the world-line integral, Eq. ( |T6l ). describes nothing but the 
proper-time evolution of a charged particle in the effective Abelian field H, governed by the 
Lorentz force oc (z x H). The world-line integration, which is equivalent to the "summation" 

n 

over all Landau levels, leads to the bosonic Euler-Heisenberg Lagrangian [19|, namely 

1 Hs 



I 



(19) 



(47rs)3/2 smh{Hs)' 

The /i-integration can be performed analytically. We further set Nc = 3 and denote 

e = Cg', 

suppressing for brevity the temperature dependence of g{T). Rescaling h = H/^'^^ we obtain 

-i4 ^ roo 



"^^smhiehs) Jo '^^(16A2/i2 + 1)2- 



(20) 
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To the leading order of the ^-expansion, sinh(^^/is) ~ C'^hs, and one can use the value of the 
integral 

dhh^ / d\- 



(16A2/i2 + 1)2 256 
to recover the Stefan-Boltzmann result: 

n=l •'^ 

The next term of the ,^-expansion of Eq. (l20l) should correspond to the g^-i^Yva in Eq. (fT4l) . 
To derive it, we should approximate ^xr^i^hs) ~ (^hs\l + (,^^/;,s)^/6], that yields 



^ „=1 ^0 ^ Jo Jo 



The /i-integration in this formula can be performed analytically, so that we obtain 
Fp T ~ —y / — e"^ / dX —. 

To the order C(^°), this formula recovers again the Stefan-Boltzmann free-energy density, 
Eq. (tai). The 0(^^)-term of interest can be obtained by approximating the sum over winding 
modes by the first two terms: 

ds f , r ^^ A'^'e-^ 17^ 277r3/2 ^ ^ 



— e-s+e-^ / rfA ^ = ^^ r77-C + C(r)- 

This yields the desired free-energy density 

F^^T)^-T'(^^-^■Cc/(T)y (22) 

which can be compared with the known result, Eq. (fT4l l. Comparing the 0{g^)-tevms, 
^ ~ 1.75 and ^ ~ 1.72, we conclude that the approximation of the full sum over winding 
modes by the (ra = 1)- and the {n = 2)-terms is indeed very good. Comparing the 0{g^)- 
terms, we finally obtain: 

We set the function C{T) to this constant value at the highest temperature of 4.54Tc achiev- 
able by the lattice simulations. With the decrease of temperature, perturbation theory 
becomes less and less applicable. We impose this fade of perturbation theory by the gradual 
nullification of the function C{T). The full free-energy density, corresponding to the spatial 
Wilson loop (fTTI) . will be calculated in the next Subsection. 
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B. The full free-energy density 



To calculate the world-line integral with the full Wilson loop, Eq. (fTTil . we first transform 
in the exponent to S by introducing an auxiliary integration as 



{W[z])=e 



dX 

7^ 



exp 



r 

4A 



+ a]S 



where 



7 = C(T)/(T)T, 

a = (j{T). By introducing one more such integration, we further transform 5* to S*^: 



{W[z]) 



vrA Jo 



exp 



We can now apply the Hubbard-Stratonovich trick, Eq. pSjl . and integrate over /z analyti- 
cally. With the use of Eq. (fT3l) . this yields 

„2 



{W[z]) 



TT 



5/2 



dX 

7=x 



d-^H- 



,iHf 



H2 + 



(24) 



The world-line integration can again be performed by using Eq. ([T9ll . that finally yields the 
full free-energy density 



X 



n=l 



F{T) = 
l^"" " X 7XV4A 



vr7/2 
+ a I e 



X 



r^ff gVsinh(gg) 







2 • 



(25) 



The remaining ordinary integrations over H, A, and s will be performed numerically in the 
next Section. 



III. NUMERICAL EVALUATION 

We fix now the temperature dependences of various quantities entering the calculation. 
We adopt the perturbative two-loop running coupling at finite temperature [3|] 



. 2.0 1,4 . I in (2 ,„ I) . whe.e . ii^, . f (^V , ,26) 



and Nc = 3 for the case under study. Furthermore, we assume the value A = 0.104Tc 



a. 
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At temperatures smaller than the temperature of dimensional reduction, T < T^, all 
dimensionful quantities can be approximated by their zero-temperature values We use 
the value = 2Tc, as suggested by the same lattice data sl- The temperature-dependent 



string tension in the adjoint representation of interest reads 



(^(T) = -(To 



1 at T^<T < T,, 



at T > 



(27) 



where o"o = (440 MeV)^ is the zero-temperature value of the string tension in the fundamental 
representation. The coefficient 9/4 stems from the Casimir-scaling hypothesis for the string 



25| to be a good 



tensions in various representations. Casimir scaling has proven in Ref. 
approximation also for the averaged Polyakov loop in the adjoint representation, (L(T)). 
That is, this average can be approximated as {L{T)) ~ (Lf(T))^''^, where (Lf(T)) is the 
averaged Polyakov loop in the fundamental representation, which has been simulated on the 
lattice in Ref. 8|. Since the (A^^ = 4)- and (A^^ = 8)-data from Ref. Q| follow the same 
pattern, we fit them both by the four-parameter function, which we use in the subsequent 
calculation: 



(Lf(r)) = 1.11 



1.53 







1 + exp 


1.92 - l) 



(28) 



In Fig. [2], we plot both the lattice data on the averaged Polyakov loop (Lf(T)) and the 
results of the fit. For comparison, in the same Fig. O we present the results of a simpler fit 
for the (A^^ = 8)-data alone [231 

1 



(Lf(T)) = exp -- [-0.23 + 1.72(T,/T)2] 



(29) 



and observe a bit less good agreement. Note that, at T ^ Tc, our fit for (Lf(T)) and the fit of 



Ref. 



27|| approach unity from above, in accordance with the positivity of the C((?^) -correction 



in Eq. 



22 



2^. Furthermore, to estimate the effect of the second winding mode, we 
use for the averaged square of the Polyakov loop in the adjoint representation the following 
approximation: (L^(T)) ^ (L(T))^. Note that the freedom in choice of the renormalization 
scheme results in the multiplication of the renormalized Polyakov loop by a temperature- 
dependent factor. Such a multiplication is equivalent to the adding of some constant A to 
the static quark-antiquark potential. We use here the values of the renormalized Polyakov 
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loop in the fundamental representation from Ref. [8|. In that paper, A is fixed to by 
demanding the static potential at T = to have the form V{r) = dor — 7r/(12r) at r > ro, 
where the value of cxo is given right after Eq. (l27l ). and ro — 0.5 fm is the so-called Sommer 



scale 



28l |. Finally, to minimize the number of free parameters, we assume that the function 



C{T) rises linearly from to the derived value (i23ll when T grows from T^, to 4.54Tc, that 



is the maximal temperature at which the lattice values 
Explicitly, the assumed function C(T) has the form 



3| for the pressure are available. 



at T, < T < T,, 

0.24 ■ j^^^fzr: at n<T < 4.54T,. 

Using the above parametrizations, we calculate numerically the pressure of valence gluons 

p.,i(r) = -F(r), 



with their full free-energy density given by Eq. (l25l) . This further yields their contribution 
to the interaction measure 

{e - 3p)vai = T-g^ - 4pvai. 

Additionally, we take into account stochastic background fields, which provide spatial con- 
finement of valence gluons in our model. Their contribution (e — 3p)stoch to the interaction 
measure is evaluated in Appendix A. We find that {e — 3p)stoch is nonvanishing and nega- 
tive definite at temperatures T > T^,, where its absolute value amounts to a few percent of 

{e - 3p)val- 

In Fig. [31 we plot the full interaction measure 

g - 3p _ {e - 3p)vai + (g - 3p)stoch 

and compare it with the lattice results from Ref. Q]. In Fig. [H we do the same for the full 
pressure, 

P_ _ Pval + Pstoch 
2^4 ~ 2^4 

Given all the approximations adopted, one can conclude that, for both the interaction mea- 
sure and the pressure, the agreements with the corresponding lattice results are rather good. 

Finally, one can define separately the perturbative and the nonperturbative contributions 
to the pressure. The perturbative contribution is given by Eq. ( l22i l with the constant C 
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Figure 2: The averaged Polyakov loop (Lf(r)) according to Ref. 
the fitting curves according to Eqs. (f28]l and ([29]) . 



(courtesy of P. Petreczky) and 



replaced by the full function C(T) given by Eq. ([30l) . Explicitly, this contribution reads 



Ppeit 



17 



27 



7r2 Stt ■ 6V4 



C(T)/(r). 



The nonperturbative contribution is accordingly defined as a difference 

Pnonpert P Ppert 



2^4 



rp4 



It includes the pressure produced by stochastic background fields and the nonperturbative 
part of the pressure produced by valence gluons. In Fig. [5l we plot Ppert/T^, Pnonpen/T'^, 
and once again for comparison the lattice data on p/T^. 



IV. SUMMARY AND CONCLUDING REMARKS 

The aim of the present paper has been an analytic calculation of the pressure p{T) and 
the interaction measure {e — 3p)/T'^ of the gluon plasma in SU(3) YM theory. Our main 
phenomenological input is confinement of valence gluons at large spatial separations and 
their perturbative interaction at small separations. The corresponding form of a spatial 
Wilson loop, Eq. (fTTI) . accounts for these interactions. 
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Figure 5: Separately shown are the nonperturbative and the perturbative contributions to the ratio 
p/T^. The lattice data, as in Fig. HI are shown once again for comparison. 

One can prove numerically that, at temperatures T > ST^, the lattice-simulated pressure 
of the gluon plasma is by more than 95% described by the Stefan-Boltzmann law along with 
the perturbative interaction between valence gluons, Eq. ( fT4l l. This fact enables us to fix 
the limiting high-temperature value C of the function C(T), which is used to parametrize 
the perturbative interaction in terms of spatial Wilson loops. The value of C, Eq. (i23l l. 
is obtained by analytically calculating the world-line integral with such a Wilson loop and 
comparing the obtained free-energy density with the known result (fT4l) . 

At high temperatures, T > 2Tc, valence gluons in our model interact both linearly at 
large spatial separations and Coulomb-like at small separations, while at T < 2Tc only the 
large-distance linear part remains. With these interactions taken into account, we manage 
to analytically calculate the world-line integral, and obtain the contribution of spatial dy- 
namics to the pressure at any temperature, Eq. (l25ll . Additionally, we account also for the 
interaction of valence gluons with the 744-gluons via the Polyakov loop, which is taken as a 



fit to the lattice data 
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25|. The critical behavior of the Polyakov loop is crucial for the 



nullification of the pressure at T = Tc and a peak of the interaction measure near Tc (cf. 
Refs. 0, Q)- 

Furthermore, stochastic background fields, that can be responsible for the 
spatial confinement of valence gluons, produce their own contributions to the interaction 
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measure and the pressure, which are also taken into account [cf. Eqs. (l33ll and (l34ll ]. In 
Figs. [3] and [U we compare our results for these two thermodynamic quantities of the gluon 
plasma with the corresponding lattice data. The observed reasonable agreement suggests 
that spatially confined valence gluons, along with the stochastic background fields, could 
indeed be the right degrees of freedom for the description of the gluon-plasma thermody- 
namics. 

Last but not least, we mention that, on the technical side, our paper advanced to an- 
alytically calculate the world-line integral with the minimal-area surface for each trajec- 
tory involved in the integration. Of course, it became possible only by virtue of a certain 
parametrization of such a surface in terms of the corresponding trajectory. Namely, we first 
assume that the surface has a parasol-shaped form, so that its area is 5 ^ | dT\z x z|, and 
use the approximation of non-backtracking planar trajectories 5* ~ | \Jq dT{z x z)|. With 
the Hubbard-Stratonovich transformation (fTSll . furthermore it is possible to single out the 
trajectory-dependence in the form exp[|H dr^z x z)] , where H is an auxiliary constant 
Abelian magnetic field. For such a field, the world-line integral becomes the known Euler- 
Heisenberg Lagrangian, Eq. f fTOD . It looks promising to generalize this parametrization to 
the 4D case, that amounts to additionally introducing an effective Abelian electric field. 
Such an approach would enable one to systematically account for effects of confinement in 
various Feynman diagrams (cf. Ref. |^). Work in this direction is in progress. 

Finally, let us mention other recent studies, where the same thermodynamic quantities 
have been calculated either by analytic methods, perturbatively in js^ and nonperturba- 



tively in 



3l| ; or on the lattice [32|, l33||. In particular, Ref. [33|| generalizing the interaction 
measure obtained in Ref. to A^'c > 3 states that e — 3p is approximately independent of 
for Nc < 6. This observation, once receiving confirmation by further lattice measurements, 
deserves theoretical interpretation. 
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Appendix A. Contributions of stochastic background fields to the pressure and the 

interaction measure 

The contribution of stochastic background fields to the interaction measure can be cal- 
culated directly, by the known formula, which assumes no valence gluons js^ 

{e - 3p)„ovai.gi. = ^ [{i9H-r)r " {i9m')o\ (31) 

with b = + y in the two-loop approximation adopted [cf. Eq. (1261) ]. and Nc = 3. 
The contribution of the chromo-magnetic condensate to the trace anomaly in the vacuum is 
subtracted. It amounts to a half of the full contribution, while the other half is represented 
by the chromo-electric gluon condensate {{gEiy)rp, whose "evaporation" at T = Tc leads 
to deconfinement {ill, Q, 21\- Thus, {{gH^Y)^ = | ((fi'-^/^j.)^)o, where the value of the 



full condensate consistent with the known values of the string tension and of the vacuum 
»n*«o„ length fl . ^ 3.55 GeV QQ. The te.pe.a..e dependence of 

the chromo-magnetic condensate can be modeled by a formula similar to Eq. ( 1271 ). 



1 at Tc < T < T,, 



at T > 



According to Eqs. (ISTll and (l32l) . stochastic background fields alone produce no contribution 
to the interaction measure at Tc < T < T^, due to the constancy of {{gH^ Y)j, at these 



temperatures [39||. At T > T^, this contribution starts appearing, but at the highest tem- 
perature under consideration, T = 4.54Tc, it must vanish again, as we assume that at this 
temperature all thermodynamic quantities are completely saturated by perturbative inter- 
actions of valence gluons. Since the increase of the perturbative contribution is given by 
the function (l30l) . the simultaneous decrease of the nonperturbative contribution should be 
described by the function [C-C{T)]/C, where C = C(4.54re) ~ 0.24, Eq. ((23]). Hence, the 
contribution of stochastic background fields to the interaction measure can be parametrized 
as 

(e - 3p)stoch = ■ [{{9Ht)\ - {{gH^)\] . (33) 
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The contribution to the pressure produced by stochastic background fields, Pstoch, can be 
obtained by integrating the formula 

stoch 

from T = T^:, below which {e — 3p)stoch = 0, on. This contribution thus reads 

Pstoch(T) _ /•^^^, (£-3p)stoch(T0 ^^^^ 

As follows from Eqs. (I32l) -(l34ll. both {e — 3p)stoch and Pstoch are negative definite, leading to 
the decreases of the full {e — 3p) and p by a few percent at T^.. < T < 4.54Tc. 
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